International  Journal  of  Electronics  and 
Communication  Engineering  (IJECE) 
ISSN(P):  2278-9901;  ISSN(E):  2278-991X 
Vol.  6,  Issue  4,  Jun  - Jul  2017;  1-8 
© IASET 


IASET 


International  Academy  of  Science, 
Engineering  and  Technology 

Connecting  Researchers;  Nurturing  Innovations 


RATE  ANALYSIS  OF  MASSIVE  MIMO  SYSTEM  USING 
STOCHASTIC  GEOMETRIC  MODEL 

SHICHAOFEN G,  SIHANLUO  & TIANLICUI 

Beijing  University  of  Technology,  Beijing,  P.  R.  China,  China 


ABSTRACT 

The  rate  performance  of  massive  multiple  input  multiple  output  (MIMO)  system  has  been  analyzed  in  the 
hexagonal  model  with  deterministically  placed  users,  which  is  too  idealized  and  lacks  rigorous  theoretical  analysis. 
This  paper  focuses  on  the  downlink  rate  performance  analysis  on  the  massive  MIMO  system.  Two  general  stochastic 
geometric  models  are  employed  owing  to  their  tractability,  i.e.,  Gilbert  disk  model  and  Voronoi  tessellation  model. 
Expressions  for  downlink  average  ergodic  rate  are  derived  in  the  general  case  and  then  simplified  in  the  special  case  with 
definite  path  loss  exponent.  The  simulation  results  of  both  models  show  that  increasing  the  number  of  available  pilots  in 
each  cell  and  the  ratio  of  base  stations  (BS)  density  to  user  equipment  (UE)  density  can  improve  the  downlink  average  rate 
per  UE.  Furthermore,  we  compare  the  downlink  rate  results  in  these  two  models  and  observe  that,  pilot  contamination  in 
the  Gilbert  disk  model  is  severer  than  that  in  the  Voronoi  tessellation  model. 

KEYWORDS:  Massive  MIMO,  Stochastic  Geometric  Model,  Pilot  Contamination,  Average  Achievable  Rate,  Wireless 
Network 

I.  INTRODUCTION 

In  massive  multiple  input  multiple  output  (MIMO)  system,  the  base  station  (BS)  employs  antenna  arrays  with  an 
order  of  magnitude  more  elements  than  the  total  number  of  serving  single-antenna  user  equipments  (UEs).  In  recent  years, 
massive  MIMO  attracts  broad  attention  because  it  can  improve  the  transmission  rate  significantly  and  save,  transmit  power 
by  an  order  of  magnitude  [1]  -[3].  Massive  MIMO  has  also  been  regarded  as  one  of  the  most  promising  technologies  in 
future  wireless  communication  system. 

Average  rate  indicates  the  performance  of  wireless  communication  system.  In  massive  MIMO  system,  the  large 
scale  antenna  array  at  BS  is  assumed  to  be  equipped  with  unlimited  numbers  of  antenna  elements  theoretically  [3]  -[5]. 
Under  such  assumptions,  the  effect  of  uncorrelated  noise  vanishes  because  of  the  law  of  large  numbers.  The  remaining 
impairment  is  inter-cellular  interference  caused  by  reuse  of  pilot  sequences  in  other  cells,  which  is  called  pilot 
contamination  [6]  -[8].  In  addition,  with  unlimited  numbers  of  antennas,  the  effect  of  small  scale  fading  also  disappears. 
Accordingly,  only  the  effect  of  large  scale  fading  remains  in  propagation.  To  sum  up,  the  rate  performance  of  the  massive 
MIMO  system  depends  on  pilot  contamination  and  large  scale  fading. 

The  main  contributions  of  this  paper  are  described  below.  Firstly,  we  discuss  the  downlink  average  rate  of  a 
typical  UE  in  the  massive  MIMO  system  in  Section  III,  using  a stochastic  geometric  model.  Expressions  for  downlink 
average  ergodic  rate  are  derived  in  the  general  case  and  then  they  are  simplified  in  the  special  case  with  definite  path  loss 
exponent.  These  expressions  which  involve  quickly  computable  integrals  are  essential  to  further  study  the  downlink  rate 
performance  of  the  massive  MIMO  system.  Secondly,  in  Section  IV,  our  analysis  demonstrates  that  the  downlink  rate 
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performance  per  UE  can  be  improved  by  increasing  the  ratio  of  BS  density  to  UE  density  or  the  number  of  available  pilots 
in  each  cell.  This  is  distinct  from  the  conventional  cellular  system  where  BS  density  doesn’t  impact  interference-limited 
ergodic  capacity  per  UE  in  downlink  [9].  In  addition,  the  downlink  rate  results  of  both  Gilbert  disk  model  and  Voronoi 
tessellation  model  are  compared.  The  reason  for  their  difference  is  analyzed  in  detail.  Simulation  results  show  that  the 
average  pilot  reuse  probability  in  the  Gilbert  disk  model  is  greater  than  that  in  the  Voronoi  tessellation  model. 

II.  SYSTEM  MODEL 


The  cellular  system  consists  of  BSs  arranged  according  to  the  homogeneous  Poisson  point  process  (PPP)  with  the 
density  A/,.  Each  BS  is  equipped  with  a large  scale  antenna  array,  containing  M (M  ->  oo)  antenna  elements.  The  density  of 
UEs  is  A u.  The  network  adopts  the  A-frequency  reuse  scheme  where  A = 1,2,3  etc.  And  there  are  P orthogonal  pilot 
sequences  available  in  each  cell,  which  means  that  the  maximum  number  of  UEs  served  by  one  BS  is  P.  All  BSs  share 
these  finite  pilots.  In  addition,  the  time  division  duplexing  (TDD)  mode  is  employed.  The  path  loss  exponent  is  a (a  >2)  in 
propagation  model. 


Channel  estimation  aims  at  acquiring  the  channel  vector  between  BS  and  its  associated  UE.  In  uplink  time  slot, 
a typical  UE  transmits  a pilot  sequence  to  its  associated  BS.  BS  accidentally  receives  signals  from  other  UEs  which  share 
the  same  pilot.  Assume  that  UE/,,,,  represents  the  UE  which  uses  the  /nth  pilot  in  the  coverage  area  of  Idh  BS  (denoted  as 


BSD-  Then  the  estimated  channel  vector  of  UE*m  acquired  by  ITS  a-  can  be  expressed  as 


iujfcm  = ln*m+L  J»«mI(m,/)  + n*, 

l*k 


where  \ M"‘ 811 ; the  channel  vector  between  BSa  and  UE/m,  gkime  CMxl  is  the  small  scale  fading  vector,  [Skim  and 

Rkim  represent  large  scale  fading  and  distance  between  the  BSa  and  the  UE/m  respectively  and  iia  is  additive  zero-mean  noise 
vector,  i Skim  is  assumed  to  obey  the  exponential  distribution.  Such  assumption  is  possible  when  the  interference  follows 
Rayleigh  fading.  1 (m,/)  is  used  to  indicate  whether  the  mth  pilot  is  reused  by  BS/,  i.e.,  = 1,  when  mth  pilot  is  used  in 

/th  cell,  0 otherwise.  After  channel  estimation,  the  downlink  data  transmission  starts  in  the  next  downlink  time  slot. 
Assume  that  Maximum  Ratio  Transmission  (MRT)  is  employed  at  BS  which  uses  the  conjugate  transpose  of  estimated 
channel  vector  as  the  precoding  vector.  Then  the  received  data  symbol  of  UEAmCan  be  calculated  as. 


II  ykm  .. 

= — = a)1"1. 

(M 


Eh/A:mX/m  . U'ktn 


fe)  W hfj  hll,n 

2^ 0 
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— ^ Olkm  0- 
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(1) 


Where,  p is  the  downlink  signal-to-noise  ratio  (SNR ),Vkm  = y/P  S/=i  htfcmX/mI(rnJ)  f Wkm  js  the  received 
signal  of  UEAm,  x/m=  h'//mS/,„  is  the  signal  transmitted  by  BS/,  .■?/,„  is  the  data  symbol  desired  by  UE/m  and  wi<m  is  the  additive 
zero-mean  noise,  (a)  Is  obtained  by  noting  that  the  correlation  of  different  channel  vectors  tends  to  be  zero  as  M ->  oo. 
(b)  Follows  from  the  law  of  large  numbers.  In  addition,  the  effect  of  small  scale  fading  disappears  as  M ->  oo  [3],  Then  the 
downlink  signal-to-interference  ratio  (SIR)  of  UEAmCan  be  expressed  as 


SIBkm  = 


PkkmRfckm 


I km. 


OO 

I km  = 'y  ^ /3lkmRikm?'(rn’  0 
l=l,l^k 


(2) 


Impact  Factor  (JCC):  4.9564 


NAAS  Rating  3.80 


Rate  Analysis  of  Massive  MIMO  System  using  Stochastic  Geometric  Model 


3 


III.  AVERAGE  ACHIEVABLE  RATE 


In  this  section,  we  will  discuss  the  downlink  average  achievable  rate.  At  first,  we  focus  on  the  probability  of  using 
the  same  pilot  between  two  BSs  in  the  system.  The  overhead  is  not  taken  into  consideration  in  this  paper.  The  average  pilot 

reuse  probability  can  be  expressed  as7?  — ® [2T(m,  l )]  — Xm=o  [N  — n]  + En=p  ^ ^ ~ n\  where  N is  the  number  of 
UEs  in  the  coverage  area  of  a typical  BS  and  P [IV  = n]  is  the  probability  density  function  (PDF)  of  the  number  of  UEs  in 
the  coverage  area  of  a BS.  Combining  (2)  with  (2),  the  average  ergodic  rate  of  a typical  mobile  user  in  the  downlink  can  be 
calculated  as. 


-(A,,  a)  = E[ln(i  - SI  If  km)]  = 111  (l- 

= ^ E,t ftt„.  [in  (l  + l1'- 

- ( fit i'}  j £/*„,.  [(«'  l)r“]  dfcir, 

•/.'  >C  J 


(3) 


Where, (r)  is  film’s  PDF,  4'  is  the  set  of  distances  between  UE  and  its  associated  BS  and  L (•]  returns  the 
Laplace  transform.  Assume  that  the  large  scale  fading  [3kkm  follows  the  exponential  distribution  with  mean  1.  (a)  Follows 
from  the  proof  of  theorem  3 in  Appendix  C of  [9].  Rkkm’s  PDF  and  set  4'  depend  on  the  employed  stochastic  geometric 
model.  Next  we’ll  present  the  expressions  of  average  achievable  rate  using  two  different  stochastic  geometric  models. 


A.  Gilbert  Disk  Model 


Gilbert  disk  model  comes  from  Gilbert’s  random  disk  graph  [10],  which  is  a special  case  of  Boolean  models  [1 1 ]. 
In  the  Gilbert  disk  model,  the  coverage  area  of  BS  is  assumed  to  be  a disk  of  radius  R centered  at  the  BS,  which  is  shown  in 
Figure  1 . UEs  are  uniformly  distributed  on  the  circle.  As  is  shown  in  Figure  1 , such  cell  deployment  inevitably  results  in 
overlap  of  different  cells. 


But,  the  Gilbert  disk  model  ensures  that  each  BS  is  at  the  center  of  each  cell.  And  the  real  cell  coverage  area 
approximates  a circle.  The  number  of  UEs  served  by  a given  BS  is  a Poisson  random  variable  with  mean  XunR2  where  Auis 
if,)  rvr  1 (A„ttP  ) ^ A..7rf?^ 

the  density  of  UE,  i.e.,  Pv  ~ n\  ~ n\  e , where  n = 0,1,2  etc.  The  average  pilot  reuse  probability  in  this 

model  can  be  expressed  as  RD  ~ Xm=o  ~p^D  — nl  + En=p  ^ — n\  Since  UEs  are  uniformly  distributed  in  each 

cell,  the  PDF  of  Rkkmis 


Figure  1:  Top  Arrangement  of  BSs  Using  Gilbert  Disk  Model,  Each  Cell  Coverage  Area  is  a Circle 
Button:  Arrangement  of  BSs  Using  Voronoi  Tessellation  Model,  the  Cell  Boundaries  form  a Voronoi  Tessellation 
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f- Rkkm  (r)  ~ R2  ' 


0 <r<R 


'Rkkm.VJ  R2'  u — ' — , where  R is  the  coverage  radius.  Assume  that  the  variable  fhkm  follows  an 

exponential  distribution  with  mean  1,  i.e.,  fhkm~  exp  (1).  According  to  the  definition  of  the  Laplace  transform, 


£!LW-E/.»  ['■  •'*"] 


/**•[- I—) 


\>ln  J]  *'*?  (-.s.j/-?,.-;,) 

IC‘  t'\Jr 

*)<ln  -’-xp  I (^  + 1)i’(-: 

H-m- 


2\ 
— I 

V 


(4) 


Where,  <f>  denotes  the  set  of  all  BSs.  (a)  is  evaluated  by  qn=  E[I(m,/)]  and  the  substitution  /?//,„,—>  f>.  (b)  follows 
from  the  probability  generating  functional  (PGFL)  of  PPP  [12],  (c)  Is  evaluated  by  using  the  change  of  variable 

s a [3  '« v2  — > t Note  that^Bff  — ^Vd/A  represents  the  effective  interfering  BS  density.  The  interfering  BS  transmits 
signal  to  its  associated  UE  which  shares  the  same  pilot  with  the  desired  UE  in  target  cell,  (d)  Follows  from  the  definition  of 

gamma  function,  (e)  Follows  from  Euler’s  reflection  formula  ^ — sin(rz)  . Using  the  substitution 


— — r ~^Tn  C 9 (s)  = exp  (— 

sincl  « ) , the  expression  can  be  rewritten  as  km  \ J. 

rate  can  be  calculated  as 


Plugging  it  in  (3),  the  average  achievable 


TD{\b,a)  = { yT  exp  [-TD  (e‘  — l)  “ r2]  ^dtdr  = ^ f 


1 - exp  \-Td  (e*  - 1)“  It2 


(a)  a 


f 


TdR. 2 

B.  Voronoi  Tessellation  Model 


L 


j p—TuR' 


(e*  - 1)< 


dt 


xn  + 1 


-dx. 


(5) 


The  Voronni  tessellation  model  is  developed  for  multi-cell  system  in  [9].  In  the  Voronoi  tessellation  model,  each 
UE  is  assumed  to  be  distributed  according  to  an  independent  PPP  with  density  Xu  and  associated  with  its  closest  BS.  The 
cell  boundaries  form  a Voronoi  tessellation,  which  is  shown  in  Figure  1.  No  overlap  between  coverage  areas  of  different 
BSs  exists.  In  reality,  each  UE  is  associated  with  its  closest  BS  so  that  the  path  loss  between  the  UE  and  its  associated  BS 
is  the  least.  Nevertheless,  BSs  in  the  model  may  be  located  so  closely  together,  but  with  a significant  coverage  area  in 
Voronoi  tessellation  model  [9].  In  addition,  the  UE  has  access  to  real  network  if  its  received  power  is  larger  than  a 
threshold.  But  UE  in  the  model  may  be  so  isolated  that  the  received  power  from  its  closest  BS  cannot  satisfy  the 
requirement. 


According  to  [13],  the  probability  density  function  of  the  number  of  UEs  in  the  coverage  area  of  a BS  is 


Fv  [iV  = n] 


3.53-5r  (n  + 3.5)  ( Xu/Xb)n 
r (3.5)  n!  (A„/Aft  + 3.5)n+3-5; 


(6) 


where  n = 0,1,2  etc.  Then  the  average  pilot  reuse  probability  can  be  rewritten  as 
P-l  oo 

Rv  = Y,  uKV  [N  = n]  + ^ Fv  [N  = n] 

n=0  n=P  . (7) 

Since  each  UE  associate  with  its  closest  BS,  the  PDF  of  Rkkm  can  be  calculated  as  [9] 
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(8) 


Assume  that  the  variable  follows  an  exponential  distribution  with  mean  1,  i.e.,  ~ exp  (1).  Then  the 

Laplace  transform  of  Ikm  can  be  rewritten  as. 


Gj  = 


= E,;. 

/ c«t\* 

?c«iV 

— E,i, 


n* I ) 

1 + X :R~" 
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-vxp[~ l 


Where,  A vefj=  Xbijv/h  is  the  effective  interfering  BS  density.  The  technical  tools  are  similar  to  Subsection  A,  so  the 
discussion  will  be  more  concise.  According  to  Appendix  C in  [9], 


£/L((e*  - iK)  = exP  ^A<#rV  - !)“  ^ 


-1)~Z  1 + Xa/2 

The  average  achievable  rate  can  be  calculated  as, 

^(e‘-l)"/(r‘-.ri7^7Td'+l]c 


dx 


(9) 


n(Ab 


3C  roc  — Aurr 
C 


d(7rA&r*’)<lf 


= r A 

Jo  r/v-(e‘  - 1)2/“  y^prTrd-r  + A 

IV.  NUMERICAL  ANALYSIS 


df. 


(10) 


In  this  section,  the  average  pilot  reuse  probability  will  be  discussed  at  first.  Then  the  downlink  rate  performance 
of  the  massive  MIMO  system  will  be  analyzed  in  detail.  Considering  the  expressions  of  average  achievable  rate,  both  of 
them  seem  to  be  irrelevant  to  BS  density.  As  is  derived  in  [9],  the  expression  of  interference-limited  ergodic  capacity  per 
user  doesn’t  depend  on  BS  density.  However,  the  BS  density  actually  impacts  the  average  achievable  rate  by  affecting  the 
average  pilot  reuse  probability  rjo  and  t]  y i n the  massive  MIMO  system.  In  order  to  express  the  rate  in  bps/Hz,  the  results  of 
(??)  and  (??)  divide  by  ln2.  The  downlink  average  achievable  rates  in  both  models  are  plotted  with  respect  to  the  ratio  of 
UE  density  to  BS  density  A„/A/,  in  Figure  2.  Firstly,  the  average  achievable  rates  decline  with  A„/A/,.  Crowded  UEs  in  the 
coverage  area  of  each  BS  will  result  in  high  pilot  reuse  probability  which  then  leads  to  intense  interference  from  pilot 
contamination.  And  increasing  available  pilots  per  cell  can  relieve  the  pilot  contamination.  Moreover,  with  a finite  number 
of  pilots,  the  rate  performance  in  a Voronoi  tessellation  model  is  worse  than  that  in  the  Gilbert  disk  model  if  Au/A/,  is 
smaller  than  the  threshold  (approximately  equal  to  P/2).  It  means  that  the  system  in  Gilbert  disk  model  achieves  better 
performance  if  dense  cells  are  deployed.  But  for  Au/A b larger  than  that  threshold,  the  rate  performance  in  a Voronoi 
tessellation  model  exceeds  that  in  the  Gilbert  disk  model.  That’s  because  the  system  in  Gilbert  disk  model  leads  to  higher 
average  pilot  reuse  probability  if  A Ll/Xb  is  larger  than  P/2.  In  the  case  A = 1,  The  lower  limit  of  downlink  average  achievable 
rate  in  the  Voronoi  tessellation  model  is  2.15  bps/Hz,  while  that  in  the  Gilbert  disk  model  is  2.02  bps/Hz. 
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?wAb 

Figure  2:  Average  Achievable  Rate  r D and  tF  vs.  Xu/Xb  for  the  number  of  Available  Pilots,  P =10,  30 
Respectively,  where  the  Frequency  Reuse  Factor  A=1  and  Path  Loss  Exponent  a =4 

V.  CONCLUSIONS 

This  paper  discusses  the  downlink  rate  performance  of  the  massive  MIMO  system.  Two  general  stochastic 
geometric  models  are  employed,  i.e.,  Gilbert  disk  model  and  Voronoi  tessellation  model.  Expressions  of  downlink  average 
ergodic  rate  are  derived  in  these  two  models,  respectively.  Simulation  results  show  that  the  downlink  average  rate  in  the 
massive  MIMO  system  increases  with  the  ratio  of  BS  density  to  UE  density  and  the  number  of  available  pilots  in  each  cell. 
In  addition,  the  rate  performances  of  these  two  models  are  compared.  Simulation  results  also  show  that  the  average  pilot 
reuse  probability  in  the  Gilbert  disk  model  is  greater  than  that  in  the  Voronoi  tessellation  model. 
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